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FIN-S is a SUPG finite element code for flow problems under active 
development at NASA Lyndon B. Johnson Space Center and within 
PECOS 

► The code is built on top of the libMesh parallel, adaptive finite element 
library 

► The initial implementation of the code targeted supersonic/hypersonic 
laminar calorically perfect gas flows & conjugate heat transfer 

► Initial extension to thermochemical nonequilibrium about 9 months ago 

► The technologies in FIN-S have been enhanced through a strongly 
collaborative research effort with Sandia National Labs 

NASA has allowed me to work here with the PECOS team since 
September 

FIN-S background and high-level overview was first presented to the 
DOE review team in October 

This talk will highlight some of new capabilities and discuss ongoing 
efforts 
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^ Software Engineering 

& Physical Modeling 

® Governing Equations 
• Thermochemistry 
® Turbulence Modeling 

Q Results 

® Viscous Reacting Flow 
® Adaptive Mesh Refinement 
® Turbulent Flow 

O Related Efforts & Ongoing Work 

® High-Temperature Thermochemistry 
® Verification 
® Near-term Effort 
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Software Engineering 


Development Environment 

• Integration into PECOS Redmine development environment 

► Source tree now housed under PECOS svn repository 

► Redmine ticket system is being used to track feature requests, 
bugfixes, etc. . . 

► Automatic Buildbot regression testing 

• Doxygen-based source code documentation 

• Rigorous modeling document 

• Example suite, unit tests, regression tests 

• GNU automake build system 
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FIN-S Code Reuse and Dependencies 

• autoconf 

• automake 

• libtool 

• Boost 

• Cantera 


• libMesh 
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Physical Modeling 


Governing Equations 


Governing Equations 

• Extension from a single-species calorically perfect gas to a reacting 
mixture of thermally perfect gases requires species conservation 
equations and additional energy transport mechanisms 

%+V-(p u) =0 
^ + V • 0 ouu) = -VP + V • r 

^ + V • (pHu) = - V ■ q + V • (tv) 
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Physical Modeling 


Governing Equations 

Governing Equations 

• Extension from a single-species calorically perfect gas to a reacting 


mixture of thermally perfect gases requires species conservation 
equations and additional energy transport mechanisms 


dp. 

dt 

dpu 

~dt 

dpE 

~df 


-Q^r + V ■ (p a u) = V • (p'Dg'Vcg) + uj s 
+ V • ( puu ) = — VP + V • t 

( ns 

p'^2 / h s V s Vc s 

s = 1 
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Physical Modeling 


Governing Equations 

Governing Equations 

• Extension from a single-species calorically perfect gas to a reacting 


mixture of thermally perfect gases requires species conservation 
equations and additional energy transport mechanisms 

^ + V ■ (p a u) = V • ( pV s Vc s ) + uj s 

dpu 
~dt 
dpE 

~df 


+ V ■ ( puu ) = —VP + V • t 

/ ns 

+ V • (pHu) = -V • q + V • (tu) + V • ( p ^ h s V s Vc s 


S=1 


• Problem class may also require a multitemperature thermal 
nonequilibrium option 


dpey 

dt 


+ V • (peyu) = -V • q v + V • ( p ^ Vc s j + cL 




s= 1 


July 21 ,2010 7/28 


Benjamin S. Kirk 


Fully-lmplicit Navier-Stokes 




Physical Modeling 


Thermochemistry 


Thermodynamics & Transport Properties 

• Thermochemistry models have been extended for a mixture of 
vibrationally and electronically excited thermally perfect gases 

e int =e trans + e rot + e vib + e elec + ^0 

ns 

= 5>ef ns (T)+ £ c,e?(T) + 

8=1 s=mol 

ns ns 

£ c.ef (T v ) + J2 c seT C (T v ) + Y / Csh° s 

s=mol s=l s=l 

Here we have assumed that T trans = T rot = T and T vib = T elec = T v 

• The transport properties have been extended as required 

► Species viscosity given by Blottner curve fits 

► Species conductivities determined from an Eucken relation 

► Mixture transport properties computed via Wilke’s mixing rule 

► Mass diffusion currently treated by assuming constant Lewis number 
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Physical Modeling 


Thermochemistry 


Chemical Kinetics 

• We consider r general reactions of the form 


N 2 + X + 

n 2 + o^no + n 


• The reactions are of the form 


7Z r 



Par 


kf r n 

S = 1 



OL sr 


where a sr and /3 sr are the stoichiometric coefficients for reactants and products 
• The source terms are then 


nr 

UJ S — ^ ^ (c^sr ^sr) ij^br ^/r) 

r= 1 
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Physical Modeling 


Thermochemistry 


Kinetic Rates 

• The forward rate coefficients are defined with a modified Arrhenius 
law as a function of some temperature T 

kfr (f ) = CfrT* exp (-E ar /Rr) 

where the rate constants are determined empirically. 

• The corresponding backward rate coefficient can be found using the 
principle of detailed balance and the equilibrium constant K eq 


• In thermal equilibrium T = T. We are currently using Cantera in this 
regime. 

• In thermal nonequilibrium T = T(T, T v ) and typical hackery ensues. 
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Physical Modeling 


Turbulence Modeling 


Turbulence Models 

• Use standard closure assumptions and eddy viscosity models 

• Spalart-Allmaras: p, t = pv sa f v 1 

^7^+V ■ (pUVsa) = C bl S s& pU sa - C w if w p (^) 

1 C/)9 

+ -V • [(/i + pv sa ) Vu sa ] H pVv sa • Vu sa 

(j a 

• k-cj (1988): fit = pk/uj 

^ + V • (/Mtfc) = r:V«- /Tpfcw + V • [(p + a*p t ) Vk] 
^ - + V • ( puu ) = : V« - fipLU 2 + V • [(/2 + apt) Vw] 

• k-Lo (2006) and SST soon to come 
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Results 


Viscous Reacting Flow 


2D Extended Cylinder 

• Laminar flow in thermal equilibrium 

• No-slip, adiabatic, noncatalytic wall 

• Chemical nonequilibrium, 5 species air (78% N 2 , 22% 0 2 ) 

Uoo = 6, 731 m /sec 
Poo = 6.81 X 10“ 4 kg/m 3 
Too = 265 K 

• Blottner/Wilke/Eucken with constant Lewis number Le = 1.4 for 
transport properties 

• Mesh, iterative convergence 

• FIN-S/DPLR comparison 

• Weak & Strong Scaling 
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Results 


Viscous Reacting Flow 
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Results 


Viscous Reacting Flow 


Mesh Convergence 
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Relative Transient Residual, \AU/At\ 


Results 


Viscous Reacting Flow 


Iterative Convergence 
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Time Step Size, At, (seconds) 




Viscous Reacting Flow 



x (m) 


Code-to-Code Comparison - 

Stagnation Line 
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Results 


Code-to-Code Comparison - 

Flank Line 
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Viscous Reacting Flow 



T(K) 



r-Stokes 
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Results 


Viscous Reacting Flow 


Speedup 
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Results 


Adaptive Mesh Refinement 



AMR - 1 3,079 node mesh, “spot on” with uniform 1 1 5,921 node mesh 


Benjamin S. Kirk 


Fully-lmplicit Navier-Stokes 


July 21, 2010 19/28 




Results 


Turbulent Flow 


Initial Turbulent Results 

• Fully turbulent flow over a flat plate 

• k-uj turbulence model; calorically perfect iV 2 ; adiabatic wall 

• Re L » 1 x 10 6 ; « 0.2 


Boundary layer profiles at trailing edge 



Code and solution verification activities are ongoing 
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Related Efforts & Ongoing Work 


High-Temperature Thermochemistry 


HTChem 

• The high-temperature thermodynamic and transport models currently 
implemented in FIN-S are one of several possible choices, and serve 
to provide the minimum set required for algorithm development 

• It is expected that these simplified models will be invalidated for 
certain problem classes and that more complex models will be 
required 

• Similar thermochemical models are required by other areas of 
PECOS research, e.g. ablation and shock layer radiation 

• The HTChem library is being developed to consolidate efforts and 
provide a common source for requisite high-temperature 
thermochemistry and transport property data 
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Related Efforts & Ongoing Work 


Verification 


Manufactured Analytical Solution Abstraction Library 

• Dearth of exact solutions necessitates method of manufactured solutions 

• Some manufactured solutions exist for the calorically perfect Navier-Stokes 
equations 

► Developed in large part by Sandia National Labs 

► Specific solutions for field, boundary condition order-of-accuracy verification 

• Existing solutions provide a necessary but not sufficient test suite 

► Will need to develop many more solutions to verify reacting flows with complex 
transport models 

• Manufactured solutions are a valuable resource that should be accessible to 
anyone 

• PECOS is developing the Manufactured Analytical Solution Abstraction 
(MASA) library to provide well-defined manufactured solutions and source 
terms for a range of physics applications 
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Manufactured solutions are being constructed and will be incorporated into the 
FIN-S regression test suite 
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Related Efforts & Ongoing Work 


Manufactured analytical solutions (used by Roy, Smith, and Ober) for each 
one of the primitive variables in Navier-Stokes equations are: 

P {x, y)= p 0 + p x sin ^ H ) + Py cos ^ ™ j , 
u (x, y) =UQ+ u x sin ^ — — — j + u y cos ^ ^ — J , 
v (x, y) = v o + cos ^ J + % sin ^ ^ J , 

P (*, y) = Po + Px cos + Py sin 
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Related Efforts & Ongoing Work 


The method of manufactured solutions applied to Navier-Stokes equations 
requires modifying the governing equations by adding a source term to the 
right-hand side of each equation: 

9(p) djpu) d{pv) _ 

dt dx dy p 

d(pu) | d{pu 2 +p - t xx ) , d(puv-T xy ) _ 
dt + dx + dy lu 

d(pv) t d(pvu - Ty X ) 1 d(pv 2 +p- Tyy) _ ^ 

o, - ! - r\ “r r\ Wv 

ot ox ay 

d(pe t ) i d(pue t + pu - ut xx - vr xy + q x ) d(pve t + pv - ur yx - vr yy + q y ) _ 
dt dx dy ~ Vet 

so the modified set of equations has a known, analytical solution. 

Symbolic representations of requisite source terms and C-source code 
have recently been generated for 2D and 3D calorically perfect gas flows. 
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Related Efforts & Ongoing Work 


- ^-(TO H 




H) [ w (SjS) + ^ («)+*] 

- te(Mi) [ % c«(!f5) + v, S m(M) +«,]“ 

:?:^^;;: ( ( ?j::li:: ( (?i :: ( ^:.1 

• - - »(S-2) (*.»(«2) + W^H 


t[.,-(?f!) + %-(»)H’] 



_ ...... 

+ [h»{if£) +*-(*$*) H’ + 3 [—(*£=) Hi H} 
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Related Efforts & Ongoing Work 
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Related Efforts & Ongoing Work 


Near-term Effort 

Additional Focus Areas 

O Physics Modeling 

► Weakly Ionized Flows 

► Surface Catalycity 

► Additional Boundary Conditions 

0 Coupling 

► Radiation 

► Ablation 

0 Adjoints 

► Sensitivity analysis 

► Adaptivity 

O Scalability 


• Push range of applicability of code through internal NASA-JSC use 
this summer 

• Perform PECOS full-system simulations using FIN-S as part of year 3 
deliverables 
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Related Efforts & Ongoing Work 


Near-term Effort 


Thank you! 
Questions? 
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